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ABSTRACT 

The theoretical intensities of the soft X-ray Fe^^"*" hnes arising from 21-31' tran- 
sitions are reexamined using a three-ion colhsional-radiative model that includes 
the contributions to line formation of radiative recombination (RR), dielectronic 
recombination (DR), resonant excitation (RE), and inner-shell collisional ioniza- 
tion (CI), in addition to the usual contribution of collisional excitation (CE). 
These additional processes enhance mostly the 2p-3s lines and not the 2p-3d 
lines. Under coronal equilibrium conditions, in the electron temperature range of 
400 to 600 eV where the Fe^^^ line emissivities peak, the combined effect of the 
additional processes is to enhance the 2p-3s lines at 16.78, 17.05, and 17.10 A, by 
~ 25%, 30%, and 55%, respectively, compared with their traditional, single-ion 
CE values. The weak 2p-3d line at 15.45 A is also enhanced by up to 20%, while 
the other 2p-3d lines are almost unaffected. The effects of DR and RE are found 
to be dominant in this temperature range (400 - 600 eV), while that of CI is 
3% at the most, and the contribution of RR is less than 1%. At lower tempera- 
tures, where the Fe^^"*" / Fe^'^+ abundance ratio is high, the RE effect dominates. 
However, as the temperature rises and the Fe^^"*" abundance increases, the DR 
effect takes over. The newly calculated line powers can reproduce most of the 
often observed high values of the (Iait.os + Iait.io) / Iai5.oi intensity ratio. The 
importance of ionization and recombination processes to the line strengths also 
helps to explain why laboratory measurements in which CE is essentially the 
sole mechanism agree well with single-ion calculations, but do not reproduce the 
astrophysically observed ratios. 

Subject headings: atomic processes — line: formation — X-rays: general — 
techniques: spectroscopic 
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1. INTRODUCTION 

The soft X-ray emission-line spectrum of Ne-like iron (Fe^^+, Fe XVII) is one of the most 
extensively investigated spectroscopic systems in astrophysics. Owing to the high abundance 
of Fe and to the closed electronic shell structure of the Ne-like ground configuration 2s^2p^, 
Fe^^+ forms over a wide range of temperatures {kT — 100 - 1000 eV). Hence, many spectra 
of hot astrophysical objects exhibit prominent Fe^^"^ lines. A better understanding of the 
spectrum of Fe^^"^ is even more important in light of the multitude of high-resolution X-ray 
spectra obtained with the grating spectrometers on board the Chandra and XMM-Newton 
observatories. The brightest lines of Fe^^"*" in the X-rays arise from the 2p-3d and 2p-3s 
transitions, which appear roughly at 15 and 17 A, respectively. The analysis of X-ray spectra 
emitted during solar flares has shown the potential of these lines to serve as diagnostic tools, 
mainly for the diagnostics of electron temperature (Raymond & Smith 1986; Smith et al. 
1985), electron density (Schmelz, Saba, & Strong 1992; Waljeski et al. 1994; PhiUips et al. 
1996), and elemental abundances (Waljeski et al. 1994). 

However, the usefulness of the Fe^^"*" lines has been hampered by two persistent problems, 
which emerge from extensive comparisons between observations, laboratory measurements, 
and theoretical models. The first is that the strongest hne at 15.01 A arising from the 
2p^ ^So - 2p^3d ^Pi transition appears to be significantly weaker than its predicted relative 
intensity. In particular, the intensity ratio of that hne to the 2p^ ^Sq - 2p^3d ^Di line at 
15.26 A, commonly labeled I3C /ho (Parkinson 1973), is consistently lower than the ratio 
calculated for a wide range of plasma parameters. This has led to the suggestion that 
the high oscillator-strength line at 15.01 A is quenched by resonant scattering (Rugge & 
McKenzie 1985; Schmelz et al. 1992; Waljeski et al. 1994; Schmelz et al. 1997; Phillips et 
al. 1996, 1997; Bhatia & Kastner 1999; Saba et al. 1999). The continuing discrepancies 
between observations and theory have prompted several laboratory studies in which the 
plasma elemental composition, and to a large extent also the plasma conditions, can be 
controlled. In particular, the EBIT (electron beam ion trap) device offers the opportunity 
to isolate nearly entirely the Fe^^+ ions from neighboring charge states and to measure their 
atomic properties with high accuracy. The EBIT measurements of Brown et al. (1998) first 
supported the claim that the line at 15.01 A is affected by resonance scattering although to 
a lesser extent than originally believed. However, as suggested by Behar, Cottam, & Kahn 
(2001a) and later demonstrated in an EBIT measurement (Brown et al. 2001), a satellite line 
of Fe^^+ that coincides in wavelength and blends with the 15.26 A line can explain most of 
the observed /3c / ho ratios without needing to invoke resonant scattering. This conclusion 
is also supported by measurements in the PLT tokamak (Beiersdorfer et al. 2001). 

The second problem with the Fe^^"*" lines is the overall intensities of the 2p-3s lines. 
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which often appear enhanced relative to theoretical predictions, when compared with the 
intensities of the 2p-3d lines. This has led to suggestions that 2p inner-shell ionization of 
Na-like Fe^^"^ (ground configuration 2p^3s) might contribute to the population of the 2p^3s 
levels (Feldman 1995). Other explanations suggested possible contributions to these level 
populations from dielectronic recombination (DR) of F-like Fe^^+ (2p^) followed by radiative 
cascades (Saba et al. 1999; Liedahl 2000; Laming et al. 2001). An KBIT measurement that 
focused on the overall intensity ratio of the 2p-3s lines relative to the 2p-3d lines (Laming et 
al. 2001) yielded a ratio that generally agrees with theoretical predictions, but is significantly 
lower than most astrophysical measurements. This further suggested that ionization and 
recombination processes, which are by and large absent from EBIT experiments and are also 
neglected in the theoretical models, play an important role in producing the Fe^^"^ lines in 
cosmic sources. 

The recent high-resolution observations of stellar coronae with Chandra (Brinkman et 
al. 2000; Canizares et al. 2000; Behar et al. 2001a; Ayres et al. 2001; Drake et al. 2001; 
Huenemoerder, Canizares, & Schulz 2001; Ness et al. 2001; Mewe et al. 2001) and XMM- 
Newton (Brinkman et al. 2001; Giidel et al. 2001a,b; Audard, Giidel, & Mewe 2001a; Audard 
et al. 2001b) produced a variety of 2p-3s / 2p-3d intensity ratios. However, none of the 
observed ratios are fully consistent with existing theoretical models, further demonstrating 
that the above problems occur for a wide range of coronal sources. Additionally, thanks to 
the high dispersion and efficiency of the Reflection Grating Spectrometer (RGS) on board 
XMM-Newton, Fe^^+ line spectra can now be obtained not only for stellar coronae (point 
sources), but also for extended sources. In observations of supernova remnants (Behar et al. 
2001b) and elliptical galaxies (Xu et al. 2002), the relative Fe^^+ line intensities show the 
same discrepancies with available models. In fact, for the analysis of the X-ray spectrum of 
NGC 4636, Xu et al. (2002) were required to invoke the empirical Capella spectrum into their 
model, because even the most updated theoretical models could not reproduce the observed 
emission. The pervasive discrepancies between the calculated spectra and those measured 
for a wide variety of hot astrophysical sources strongly suggest that the problems lie in 
the shortcomings of the atomic models for line formation rather than in the uncertainties 
associated with the astrophysical environments. 

In an effort to improve the existing atomic models, which neglect atomic processes that 
involve neighboring ions, we construct a 3-ion coUisional-radiative model that includes levels 
of F-like, Ne-like, and Na-hke Fe. Such a model allows one to consider coUisional excitations 
as well as the effect of recombination and ionization simultaneously. The main drawback of 
a naive, all-inclusive 3-ion model, in which rate equations for the three ionization states are 
solved simultaneously, is that the fractional ion abundances, which self-consistently result 
from the model, would be clearly erroneous. In order to obtain correct results for the 
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ionization balance in a coUisional-radiative model for a wide range of temperatures, one 
would have to include more than three ionization states, as well as numerous levels for each 
charge state (to adequately account for DR), which makes this approach impractical. Since 
we are interested in the line spectrum and not the ionization balance, we can circumvent 
this difficulty by adopting the approach formerly used in Doron et al. (1998, 2000), by which 
an independent ionization balance is imposed on the level-by-level, multi-ion model. The 
details of the current 3- ion method are given in the following section (§2). In §3, we give 
the results of the 3-ion model for the line-emission of Fe^^+ and compare them with results 
obtained with a basic, single-ion model typical of those that have been used until now. 



CoUisionally ionized, astrophysical plasmas are characterized by relatively high temper- 
atures and low densities. The most common atomic model used to describe the line emission 
from a particular ionic species (say Fe^^"*") under these conditions includes only electron im- 
pact excitations from the ground state and subsequent radiative decays, either directly to 
the ground state or via cascades. In the present work, we wish to extend this standard 
model and take into account also atomic processes involving neighboring ionization states. 
Specifically, we include dielectronic and radiative recombination from F-like Fc^^+, inner- 
shell ionization from Na-like Fe^^+, and resonant excitation through doubly-excited levels of 
Fe^^"*" (i.e., dielectronic capture followed by autoionization to excited levels). We assume the 
plasma is optically thin and coUisionally ionized, i.e., photoexcitation and photoionization 
are unimportant. Assuming that the free electron energy distribution is Maxwellian, cor- 
responding to an electron temperature Te, the general set of rate equations for the density 
(population) rf-^ of an ion with charge q-\- in a level j in a steady state can be written as: 



2. 



THEORETICAL METHOD 
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where / and j" are levels pertaining to ions with charges (g — 1)+ and {q + l)+, respectively. 
The terms Q{Tf.) and Q^^{Te), respectively, represent the rate coefficients for direct and 
resonant (RE) electron-impact excitation or de-excitation. The terms A and A"" denote the 
rate coefficients for spontaneous radiative decay and for autoionization, respectively. The 
terms S{Te) stand for the rate coefficients for collisional ionization, and a^*" (Te) and a^^{Te) 
are the rate coefficients for dielectronic capture and for radiative recombination, respectively. 
The electron density is denoted by He- The entire set of equations that needs to be solved 
consists also of the appropriate rate equations for n^^~^^^ and n^^,^^^~^, consistent with the 
atomic processes that appear in eq. (1). Throughout this work, we use q — 16. Densities for 
ions other than the three considered are taken as zero. In order to obtain a unique solution, 
one of the rate equations in cq. (1) is substituted by a normalization condition, such as 
^i' ~^ ~^ ^i" ^i"^^^^ ~ ^' "^liich leaves us with N — 1 independent equations 

for the N levels included in the model. 

A straightforward solution for the set of equations (1) necessarily yields the ionization 
balance through the results for the level populations nj? n'j'^ , and n^v,^^'*^. As argued 
above, this resulting ionization balance would clearly be wrong for most temperatures, as 
the set of equations consists of only three ionization states, and practically, an insufficient 
number of levels to determine the correct ionization balance. In order to avoid this error, 
the set of rate equations is solved by imposing the following constraints: 

E-r^vE (2) 

j' i 

E</E4^'^^-A//.+i (3) 

where /g, and Z^+i represent the relevant ion fractions, and are calculated indepen- 

dently. In other words, we impose the ionization balance and only seek for the consequential 
population distribution among the individual levels. We are primarily interested in levels 
j of the ion with charge g+. In the coronal approximation, where the ions arc essentially 
in their ground state, and in the absence of mctastable levels, it is sufficient to apply the 
constraints only through the ground levels of the three charge states. Moreover, collisional 
transitions among excited levels can be neglected. For Fe^^"*", this approximation holds for 
Ue < 10^^ cm~^. For the present study, we have adopted the fractional abundances calcu- 
lated by Mazzotta et al. (1998). Table 1 presents the ion fractions of Fe^^"*" through Fe^^"*" for 
several chosen electron temperatures in the relevant range of 200 - 1000 cV. The last column 
in the table gives the combined fraction of these three ions as a function of temperature. 
These total values are provided in order to give an idea of the collective importance of these 
ions in astrophysical plasmas at each temperature. 
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For Fe-*^^"*" we include in the model the configurations 2s^2p^, 2s^2p^nZ (n = 3 to 7), 
2s2pSZ (n = 3 to 7), 2s22p^3/3/', 2s2p53/3/', 2s22p43/4d, and 2s22p^3d4Z. For Fe^^+, the 
configurations 2s^2p^ and 2s^2p"^3/ are taken into account. Finally, the ground configuration 
2s^2p®3s of Fe^^+ is also included. All of these configurations add up to 2111 levels. In fact, 
the inclusion of the RE processes (Q^^) in the set of rate equations (1) is equivalent to 
including the doubly-excited configurations of Fe^^+. Since we are not interested specifically 
in the population of the Fe^^"^ levels, we have incorporated the RE rate coefficients from 
the literature (instead of including all of the appropriate doubly-excited levels of Fe^^+), 
and by doing so, we significantly reduce the number of rate equations to be solved. The 
singly-excited configurations of Fe^^+ up to n = 7 have been shown to account for most of 
the radiative-cascades effects (Liedahl 2000) and are taken into account in all of our models. 
The Fe^^+ doubly-excited configurations that are included here explicitly produce about 50% 
of the total DR effect in the entire 200 - 800 cV temperature range. The contribution of 
higher DR resonances is obtained by means of extrapolation based on the total DR rates 
published by Chen (1988). In order to account for the radiative recombination contribution 
onto high electronic shells with n > 7, we have employed a hydrogenic extrapolation method. 

The solutions for the rate equations (1), which are the normalized level populations 
nj"*"/ '^k'n'f^, are used to obtain the line power /j^, which is defined here by: 

k 

Subsequently, the line emissivity which gives the total number of photons emitted in the 
plasma per unit time and per unit volume as a result of the radiative decay from level j to 
level i, can then be found by: 

If, = ifneApeUH (5) 

where Ap^ is the iron elemental abundance and uh is the hydrogen density. 

The basic atomic quantities used in this work were generated by means of the multi- 
configuration, relativistic HULLAC (Hebrew University Lawrence Livermore Atomic Code) 
computer package developed by Bar-Shalom, Klapisch, & Greg (2001). Resonance excitation 
rates are not calculated explicitly, but incorporated from Chen & Reed (1989) for the 2p^ to 
2p^3/ (/ = s, p) excitations of Fe^^"*". Following those authors, resonant excitation to 2p^3d 
levels is assumed to be unimportant due to their high energies and has been neglected. 
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3. RESULTS AND DISCUSSION 

We have employed the model described above to calculate the Fe^^+ level populations 
and line intensities for various temperatures in the range 200 - 1000 eV. In the following, we 
focus our presentation on the results for the seven most intense Fe^®+ lines, which fall in the 
range 13.5 - 17.5 A. The corresponding transitions, their labels, and measured wavelengths 
(A) from Brown et al. (1998) are given in Table 2. Since in optically thin plasma, the relative 
intensities of these lines show negligible sensitivity to the electron density, the present results 
are vahd for electron densities below 10^^ cm~^. 



3.1. Line Powers 

The calculated line powers [defined in §2 eq. (4)] as a function of the electron tem- 
perature for the seven transitions hsted in Table 2 are given in Figs. 1 and 2 by the solid 
curves. For comparison, we also give the line powers calculated with a single-ion model, 
i.e. a model that takes into account only direct electron impact excitations and subsequent 
radiative decays. The data points in the plots throughout this paper mark the explicitly 
calculated values, while the curves represent a spline interpolation between these points. It 
can be seen from this comparison that the 2p-3s lines (Fig. 1), as well as the weak 2p-3d 
line at 15.45 A (Fig. 2), are significantly affected by processes involving neighboring ions 
over a wide range of temperatures, while the two strong 2p-3d lines at 15.01 and 15.26 A, 
as well as the 2s-3p line at 13.83 A, are much less affected by these processes. In particular, 
the intensity of the strongest 2p-3d fine at 15.01 A is almost identical in the two types of 
models, while the intensity of the 2p-3s fine at 17.10 A roughly doubles when contributions 
from neighboring ionization states are considered. 

3.2. Contribution of the Various Mechanisms to Upper Level Population 

In order to obtain an insight into the dramatic differences between the results of the 
3-ion model and the single-ion model, we investigate the relative importance of each atomic 
process separately. This is done by running a sequence of 3-ion models, where in each run we 
turn off one of the atomic processes: resonant excitation (RE), dielectronic recombination 
(DR), radiative recombination (RR), or collisional ionization (CI). The resulting line powers 
in each of these models are presented in Table 3 for several electron temperatures in the 
range 200 - 800 eV. The results obtained with a traditional single-ion model, which includes 
direct collisional excitation (CE) alone, are also given for comparison. The first column in 
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the table indicates the electron temperature and the second column describes the model. 
The following columns give the line power for each of the seven lines. In general, it can be 
seen from Table 3 that among the additional atomic processes, which are included in the 
3-ion model, RE and DR are the most important. 

At an electron temperature of 400 eV, the power of the 17.10 A line, which is the 
most affected, receives contributions of 16%, 11%, 0.5%, and 3% from RE, DR, RR, and CI, 
respectively, which together enhance the line by 47% compared with the value obtained when 
calculating direct CE alone. At 600 eV, where the absolute Fe^^+ line powers are still very 
high, the additional processes enhance the CE value for the 17.10 A line by 57%, where RE, 
DR, RR, and CI, respectively, produce 10%, 21.5%, 0.5%, and 2.5% of the total line power. 
The other 2p-3s hues at 16.78 and 17.05 A are also enhanced by the non-CE processes, 
albeit to a somewhat lesser degree (20% - 30%). The main reason for this is DR, which 
preferentially populates the upper level of the 17.10 A hue. Additionally, CI processes to 
that level are about twice as effective as CI to any other upper level of the 21-31' transitions. 
The weak 2p-3d line at 15.45 A is also enhanced, mostly by DR (RE processes towards the 
2p^3d levels are neglected). At 600 eV, 13.5% of the 15.45 A hue power is due to DR. 

The trends of the various contributions with temperature are also interesting to note. 
The RE effect is predominant at low temperatures. It represents about a third of the 2p-3s 
line powers at 200 eV, but declines rather sharply with temperature. At 800 eV it contributes 
less than 7% to the 2p-3s line powers. The effect of DR, on the other hand, is small at low 
temperatures, but increases dramatically as the temperature increases. The enhancement 
of the 2p-3s lines due to DR reaches almost 30% at 800 eV. These two opposing effects 
are a direct consequence of the strong variations in the fractional ionic abundances with 
temperature and the fact that the RE effect is correlated with the Fe^^"*" abundance, while 
that of DR depends strictly on the Fe^^"'' abundance. As can be inferred from Table 1 
and from Mazzotta et al. (1998), the abundance of Fe^^"*" dominates at temperatures up 
to ~400 eV, while that of Fe^^"*" peaks at ~600 eV. Interestingly, through most of the 
temperature range in which the Fc^^"*" line powers are relatively high, say 300 - 700 eV (see 
Figs. 1 and 2), both RE and DR are important, which means neither can be neglected 
when accurate line powers are needed. Recent unpublished calculations by Chen & Pradhan 
(2002) claim that Chen & Reed (1989) have underestimated the RE rate coefficients. Indeed, 
higher RE rates would imply even a more pronounced enhancement to the line powers 
than demonstrated here, especially in the low temperature regime, where RE is important. 
Unfortunately, the absence of full level-by-level RE rate coefficients in that paper precludes 
us from estimating this effect quantitatively. Notwithstanding, the DR and CI enhancements 
calculated in this work are by and large independent of the RE effect. 
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It is important to distinguish between the DR enhancement of Fe^^+ hnes calculated 
in this work and the effect of Fe^^~^ DR satellite lines, which was traditionally incorporated 
in plasma models (see e.g., Raymond & Smith 1986). The DR effect, which is discussed 
here, produces the exact same parent lines of Fe^^"*" (i.e., 2p^3/ - 2p^), by means of DR 
processes that populate the upper levels of these transitions (e.g., 2p'^3/3/' - 2p^3/). This 
should not be confused with DR satellites emitted by lower charge states at slightly longer 
wavelengths (e.g., 2p^3/3/' - 2p^3Z in Fe-^^+). The latter, in many cases, can be resolved from 
the parent lines by high-resolution spectrometers. Note, however, that the radiative decays 
(2p^3Z3Z' - 2p^3Z') that populate the upper levels of the Fe^^"*" hnes, are in fact associated 
with DR sateUites of Fe^^"*". This hue emission is exphcitly included in our calculations and 
is illustrated in the spectra shown below in §3.4. 



3.3. Line Intensity Ratios 

The large impact of the additional atomic processes considered in the 3-ion model on 
some of the line intensities directly affects the line ratios and, thus, the diagnostic appli- 
cations of the Fe-^^"*" system. In Figs. 3 - 5, we present the most interesting line intensities 
normalized in most cases to the 2p-3d line at 15.01 A, which is the strongest Fe^^+ line at 
most temperatures and also the least affected by the additional processes (see Table 3). As 
in Figs. 1 and 2, the solid curves represent the results of the (complete) 3-ion model and the 
dotted curves represent the single-ion model. The most pronounced effect is found for the 
ratios involving the 2p-3s lines in Fig. 3. Since the spectral resolution of contemporary spec- 
trometers, such as the Reflection Grating Spectrometer (RGS) aboard XMM-Newton and 
the Low Energy Transmission Grating Spectrometer (LETGS) aboard Chandra, is hardly 
sufficient to resolve the two 2p-3s lines at 17.10 and 17.05 A, we also show the corresponding 
curve for the summed intensities of these two lines. The ratio involving the third 2p-3s 
line at 16.78 A is given separately in Fig. 4 to avoid overcrowding the plot. The calculated 
ratios shown in Fig. 3 shed new light on the discrepancy as yet between the observed and 
calculated values for the ratio of the 17.10 A and 17.05 A lines to the line at 15.01 A. Using 
the single-ion model, one obtains (in the relevant temperature range) a maximum value of 
1.2 for the ratio (Iait.os + Iait.io) / Iais.oi, while the observed ratios from various coronal 
sources, including many observations of the Sun, are in the range of ~ 1.5 -2.6. A list of 
various measured ratios for these hnes is presented in Table 4. The inclusion of RE, DR, 
and CI processes in the calculations by means of the 3-ion model yields values of ~ 1.4 - 
1.9 for this ratio at temperatures of 200 eV and higher, where Fe^^"*" forms predominantly. 
Even higher values, which are observed mostly in active regions on the Sun (e.g., Phillips et 
al. 1997; Saba et al. 1999) can be attained for this ratio at lower temperatures, which are 
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realistic for these environments. In the special case of spatially resolved observations of the 
Sun, resonant scattering of the 15.01 A line is also plausible. Alternatively, non-equilibrium 
ionization (NEI) conditions can also play a role in enhancing the 2p-3s / 2p-3d ratios. In 
particular, NEI conditions make the effects of ionization and recombination on line forma- 
tion more pronounced than in equilibrium and would tend to enhance the 2p-3s lines even 
further. In conclusion, the major part of the observed ratios can be explained in terms of 
the 3-ion coronal-equilibrium model. For the abnormally high ratios observed mostly in the 
Sun, resonant scattering and NEI are also plausible. 

The present results can shed more light on the intcpretation of laboratory measurements. 
Laming et al. (2001) found for the 2p-3s / 2p-3d ratios good agreement between their EBIT 
measurements and simple, single-ion models. In the EBIT experiment, the selectivity of 
the ionic charge Fe^^+ is supposed to be very high and the presence of Fe^^+ and Fe^^+ 
can therefore be minimized. In addition. Laming et al. (2001) used a relatively high-energy 
(~1 keV) electron beam. Under such conditions, the contributions of RE, DR, and CI 
are expected to be very small, and a single-ion CE model could therefore be adequate for 
reproducing the experimental results. However, in a very recent EBIT experiment, also with 
a high-energy beam (>0.8 keV), Beiersdorfer et al. (2002) measure 2p-3s / 2p-3d ratios that 
are almost a factor of 2 higher than those measured by Laming et al. (2001). These recent 
results could imply that even the basic CE models may be inadequate, or that high-lying 
resonances do contribute to RE even at these high energies, or a combination thereof. In 
any event, the effects of the additional processes, which are put forward in this work, are 
relevant to all plasmas in coUisional equilibrium and as such, they ubiquitously enhance 
the CE (-|- additional RE) rates measured with EBIT. Another indication for additional 
processes in plasma sources that are not present in EBIT experiments comes from the relative 
intensities of the 17.05 and 17.10 A lines. In the EBIT measurements of Brown et al. 
(1998) the line at 17.05 A is significantly stronger than the line at 17.10 A, while in many 
astrophysical observations the two lines have comparable intensities (e.g., Rugge & McKenzie 
1985; Huenemoerder et al. 2001; Behar et al. 2001a). This difference is consistent with our 
results as presented in Table 3, where DR enhances the 17.10 A hne much more than it does 
the 17.05 A hne. 

As can be seen in Fig. 3, for kT^ < 600 eV, the newly calculated 2p-3s / 2p-3d ratios are 
more sensitive to the electron temperature than those calculated with the single-ion model, 
while at higher temperatures, the little sensitivity that existed in the single-ion model is, 
in fact, lost in the new model. In particular, the major role played by the RE processes 
at low temperatures significantly increases the temperature sensitivity in this regime. The 
(Iai7.05 + Iai7.io) / Iai5.oi ratio decreases from 1.9 to 1.4 from 200 to 600 eV. The temperature 
sensitivity of these ratios diminishes rapidly thereafter, but in any case the absolute power of 
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these lines drops sharply beyond 600 eV (See Figs. 1 and 2). Therefore, in multi-temperature 
sources often unresolved in astrophysical observations, the high-temperature contribution to 
the integrated spectrum of Fe^^+ would be very small. In Fig. 4, one can see that the values 
obtained for the ratio Iais.oi / Iai5.26 in the 3-ion model are not very different (< 3%) than 
the single-ion model values in the entire temperature range. We note that our single-ion 
calculations for the ratio Iais.oi / Iai5.26 yield a value of 3.4. The present model includes 
configurations up to n—7. The inclusion of even higher configurations would have only a 
neghgible effect on the ratio. A value of 3.4 is about 12% higher than the value measured by 
Brown et al. (1998) in EBIT with a monoenergetic beam at 1.15 keV, and about 15% higher 
than the EBIT ratio measured by Laming et al. (2001) with a beam energy of 0.9 keV. The 
EBIT values are associated with errors of about 5%. Both measurements are just within the 
15% accuracy range expected from the CE rate coefficients produced by the HULLAC code 
for this kind of atomic system. 

The ratio of the weak 2p-3d line at 15.45 A to the 2p-3d resonance line at 15.01 A is also 
enhanced, albeit to a lesser extent than the 2p-3s lines, as can be seen in Fig. 5. The 2s-3p 
line at 13.83 A is only very slightly affected by the additional processes (Fig. 5). The 15.45 A 
line is enhanced for the most part by DR, while the effects of CI and RR are very small. 
The DR effect increases with temperature due to the increase in the abundance of Fe^^"*". 
The CI effect, on the other hand, depends both on the CI rate coefficients that increase with 
temperature, but also on the decreasing Fe^^+ abundance. Since at low temperatures, the 
rate coefficients for inner-shell ionization (CI) are very small while the Fe^^+ abundance is 
high and vice versa at high temperatures, the CI effect is supressed in ionization equilibrium 
conditions throughout the relevant temperature range. In NEI, however, the CI effect could 
be more pronounced. Note that following Chen & Reed (1989), RE processes to 2p^3d levels 
are neglected in our model. However, according to Chen & Pradhan (2002), RE processes 
do, in fact, populate the upper level of the 15.45 A line and enhance this ratio to values of 
0.08 and higher. On the other hand, the measurements of Beiersdorfer et al. (2002) show a 
value of 0.04 for this ratio, which is totally consistent with the present results. Either way, 
the DR effect on this ratio demonstrated in Fig. 5 remains. Admittedly, the ratios in Fig. 5 
could be rather hard to use for diagnostics due to their weakness in most sources. However, 
in bright sources, such as Capella for instance, the 13.83 A line could be intense enough to 
provide useful temperature diagnostics (Behar et al. 2001a). 

We point out that the present line ratios have been compared with those calculated 
by Gu (private communications) using a similar model based on the FAC atomic code (Gu 
2002). The overall 2p-3s / 2p-3d intensity ratios obtained with the two codes agree to within 
1%. The agreement between ratios for individual lines is ~5%. 



3.4. Synthetic Spectra 



As seen in §3.2 and §3.3, additional processes to direct collisional excitation (CE) have 
a major role in producing the soft X-ray lines of Fe^^+. In order to provide a more visual 
illustration of the overall effect on the spectrum, total and partial synthetic spectra calculated 
for 400 eV are presented in Fig. 6. The upper trace in the figure gives the total spectrum, 
while the lower trace shows the explicit contribution of the CE processes alone. The middle 
trace represents the complementary contribution of non-CE processes. The different plots 
in Fig. 6 are all on the same scale. The spectra clearly show the importance of the non-CE 
processes to the 2p-3s lines at ~ 17 A. The 2p-3d lines (~ 15 A) are almost unaffected 
by these mechanisms. The spectra also include lines of Fe^''"'" around 14.2 A (2p-3d) and 
16 A (2p-3s), which are produced self consistently by the model and are rather weak at 
400 eV. The line intensities for Fe^^+, however, should be viewed with caution, since the 
model for Fe^^"*" is hmited. Particularly, DR processes from Fe^^"*" are not taken into account. 
Additionally, the model does not include the DR satellite hues of Fe^^+ and therefore these 
are absent from the plots in Fig. 6. At 400 eV, these DR sateUites are expected to produce 
weak line emission mostly around 15.26 A, as mentioned above in §1. 



4. CONCLUSIONS 

We have constructed a 3- ion model to account for ionization and recombination processes 
that produce the Fe^^+ emission hnes, in addition to the coUisional and radiative processes 
among Fe^^"*" levels. It is shown that recombination (DR and RR), inner-shell ionization 
(CI), and resonant excitation (RE) processes, contribute significantly to the hne powers and 
can resolve the existing discrepancies between observations and previous calculations within 
the framework of a coronal equilibrium model. Among the additional processes considered, 
RE and DR (of Fe^''^) are found to play a major role in producing the Fe^^+ lines. The 
present calculations reproduce the systematically high values observed for the 2p-3s / 2p-3d 
intensity ratios in many astrophysical sources. Only abnormally high ratios observed mostly 
on the Sun may still require one to invoke resonant scattering of the 15.01 A line, or less 
likely transient conditions that depart from ionization equilibrium. Since non-equilibrium 
conditions are not required to explain most of the observed ratios, we defer the study of 
those cases to future work. The fact that recombination and ionization are important makes 
the present predictions dependent upon the ionization balance of Fe^^+ - Fe^''+. We expect 
similar effects of inner-shell ionization and of recombination on the 21 - 31' lines to be present 
in more highly ionized systems (Fe-'^^+ - Fe^-'^+) as well. However, a reliable assessment of 
this effect will require detailed models similar to the one employed here. The results of 
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the present, elaborate, 3-ion model should also be important for analyzing the intensities 
of other lines that arise from the decay of the 2p^3/ levels, i.e., 3/ - 3/' transitions. In the 
case of Fe^^"'", these transitions fall in the EUV range. In lighter Ne-like ions, these 3/ - 31' 
transitions are clearly observed in the UV solar spectra obtained by the Solar Ultraviolet 
Measurements of Emitted Radiation (SUMER) spectrometer on board the SOHO satellite. 

EB acknowledges ongoing work and many useful discussions with Steven Kahn on the 
complexity of the line intensity patterns of Fc^^+ and their manifestation in X-ray observa- 
tions. We arc grateful to Ming Feng Gu for useful comparisons of the present results with 
some of his own calculations. We thank Daniel Savin for reading the manuscript carefully 
prior to its submission. RD was supported by a NASA Solar Physics Guest Investigator 
Grant S137816. 
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Table 1. Ionization fractions for Fc^^+, Fe^^+, and Fc^''+ interpolated from Mazzotta et al. 

(1998) in the range 200 - 1000 eV. 



kTe (eV) Fe^5+ Fe^^+ Fe^^+ Total 



200 


0.248 


0.325 


0.0029 0.576 


400 


0.131 


0.667 


0.154 0.952 


600 


0.036 


0.353 


0.295 0.684 


800 


0.0047 


0.077 


0.137 0.219 


1000 


0.00039 


0.0088 


0.027 0.036 



Table 2. Transitions and measured wavelengths of the seven strongest lines of Fe-^^"*" in the 

range 13.5 - 17.5 A. 



Label Transition A (A) ^ 



3A 


2s22p6 


^Sq - 


2s2p63p 




13.825 


3C 


2s22p6 


^Sq - 


2s22p53d 


^Pi 


15.014 


3D 


2s22p6 


%- 


2s22p53d 




15.261 


3E 


2s22p6 


%- 


2s22p53d 


^Pi 


15.453 


3F 


2s22p6 


%- 


2s22p53s 


^Pi 


16.780 


3G 


2s22p6 


%- 


2s22p53s 


^Pi 


17.051 


M2 


2s22p« 




2s22p^3s 




17.096 



^Brown et al. (1998) 
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Table 3. Fe^^"*" line powers calculated with various 3-ion models. 



kT, (eV) 


Model 


13.83 A 
2s-3p 


15.01 A 
2p-3d 


Line power (10 
15.26 A 15.45 A 
2p-3d 2p-3d 


cm^ s ^) 
16.78 A 
2p-3s 


17.05 A 
2p-3s 


17.10 A 
2p-3s 


200 


Complete 


0.128 


2.96 


n QO'7 


U.loz 


2.39 


3.02 


2.52 




No RE 


0.128 


2.96 


U.oy ( 


U.ioz 


1.69 


2.07 


1.66 




No DR 


0.128 


2.96 


895 


1 51 


2.38 


3.00 


2.49 




No RR 


0.128 


2.96 


0.897 


0.152 


2.39 


3.02 


2.52 




No CI 


0.128 


2.96 


0.897 


0.152 


2.36 


2.98 


2.46 




CE only 


0.128 


2.96 


0.894 


0.150 


1.64 


2.00 


1.56 


400 


Complete 


2.78 


39.7 


11. 


l.oU 


25.0 


31.4 


25.8 




No RE 


2.78 


39.7 


lie 
ii.o 


i.OU 


21.3 


26.4 


21.4 




No DR 


2.76 


39.6 


11.7 


1.49 


24.2 


29.8 


22.9 




No RR 


2.78 


39.7 


11.8 


1.60 


25.0 


31.2 


25.7 




No CI 


2.77 


39.7 


11.8 


1.57 


24.6 


30.8 


25.0 




CE only 


2.75 


39.6 


11.7 


1.45 


20.1 


24.2 


17.6 


OUU 


Complete 


Old 

o.ib 


o7.o 


11.1 


1.33 


Zl.D 


ZD.O 






No RE 


3.16 


37.5 


11.1 


1.33 


19.5 


23.6 


20.3 




No DR 


3.13 


37.3 


10.9 


1.15 


20.1 


24.5 


17.6 




No RR 


3.16 


37.5 


11.1 


1.32 


21.5 


26.2 


22.4 




No CI 


3.15 


37.5 


11.1 


1.33 


21.1 


25.8 


22.0 




CE only 


3.12 


37.4 


10.9 


1.13 


17.7 


21.2 


14.4 


800 


Complete 


1.01 


10.9 


3.22 


0.366 


5.87 


7.48 


6.33 




No RE 


1.01 


10.9 


3.22 


0.366 


5.48 


6.95 


5.87 




No DR 


0.99 


10.8 


3.14 


0.292 


5.33 


6.48 


4.35 




No RR 


1.01 


10.9 


3.22 


0.361 


5.87 


7.44 


6.30 




No CI 


1.00 


10.9 


3.20 


0.365 


5.80 


735 


6.19 




CE only 


0.98 


10.8 


3.12 


0.287 


4.84 


5.78 


3.72 
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Table 4. The (Iai7.o5 +Iai7.io) /Iais.oi line intensity ratio. The reference key is: 1- Rugge & 
McKenzie (1985), 2- PhiUips et al. (1997), 3- Saba et al. (1999), 4- Canizares et al. (2000), 
5- Behar et al. (2001a), 6- Mewe et al. (2001), 7- Audard et al. (2001b), 8- Ayres et al. 
(2001), 9- Huenemoerder et al. (2001), 10- Xu et al. (2002). 



Source 


Spectrometer 


(Ia17.05 +Ia17.1o)/Ia15.0i) 


Reference 


Sun 


Solex B RAP 


1.95 


1 


Solar flares 


flat crystal 


1.93 ^ 


2 


Solar active region 


flat crystal 


2.2 ^ 


2 


Solar active region 


Hat crystal 


2.6 ^ 


3 


Capella 


Chandra HETGS 


1.67 ^ 1.65 


4,5 


Capella 


Chandra LETGS 


1.62 


6 


Capella 


XMM-Newton RGS 


1.52 


7 


HR 1099 


Chandra HETGS 


1.67 


8 


II Pegasi 


Chandra HETGS 


2.10 


9 


NGC 4636 


XMM-Newton RGS 


1.59 


10 



^averaged over many observations 
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200 300 400 500 600 700 800 900 1000 

kT^ (eV) 



Fig. 1. — Line powers of the 2p-3s transitions at 17.10 A (squares), 17.05 A (circles), and 
16.78 A (triangles). 
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200 300 400 500 600 700 800 900 1000 

kT^ (eV) 

Fig. 2. — Line powers of the 2p-3d transitions at 15.01 A (up triangles), 15.26 A (circles), 
and 15.45 A (squares), and the 2s-3p transition at 13.83 A (down triangles). 
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xn.io ;vi5.oi 

' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 



200 300 400 500 600 700 800 900 1000 

kT^ (eV) 

Fig. 3. — Intensity ratios of the two 2p-3s lines at 17.05 A and 17.10 A to the 2p-3d resonance 
hne at 15.01 A. 
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Fig. 4. — Intensity ratios of the 2p-3s line at 16.78 A and the 2p-3d hne at 15.26 A to the 
2p-3d resonance hne at 15.01 A. 
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200 300 400 500 600 700 800 900 1000 

(eV) 

Fig. 5. — Intensity ratios of the 2s-3p line at 13.83 A and the 2p-3d line at 15.45 A to the 
2p-3d resonance line at 15.01 A. 
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Fig. 6. — Synthetic Fe^^"'" and Fe^^"*" spectra calculated for an electron temperature of 400 eV. 
The upper trace is the total spectrum calculated with a 3-ion model. The bottom trace rep- 
resents the CE component alone and the middle trace includes all of the other contributions, 
namely RE, DR, RR, and CI. 



